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Abstract. - The exact elastodynamic scattering theory is constructed to describe the spectral 
properties of two- and more- cylindrical cavity systems, and compared to an elastodynamic 
generalization of the semi-classical Gutzwiller unstable periodic orbits formulas. In contrast to 
quantum mechanics, complex periodic orbits associated with the surface Rayleigh waves domi- 
nate the low-frequency spectrum, and already the two-cavity system displays chaotic features. 



Introduction. - The Gutzwiller semi-classical quantization of classically chaotic systems 
relates quantum observables such as spectral densities to sums over classical unstable periodic 
orbits [1,2]. The work presented here is a step toward a formulation of such approximate short- 
wavelength theory of wave chaos for the case of linear elastodynamics. Why elastodynamics? 
The experiments initiated in ref. [3] attain Q values as high as 5 • 10^, making spectral 
measurements in elastodynamics competitive with measurements in microwave cavities at 
liquid helium temperatures [4,5], and vastly superior to nuclear physics and room temperature 
microwave experiments for which the Q values are orders of magnitude lower, typically ~ 10^- 
10^. For elastodynamics there are only a few experimental demonstrations [6] of the existence 
of unstable periodic orbits, and no theory that would predict them. While Oxborrow et 
al. [3] measure about 10^ spectral lines, the current theory is barely adequate for computation 
of dozens of resonances. A more effective theory would find many applications such as in 
the frequency domain quality testing for small devices built from high Q materials. This 
unsatisfactory state of affairs is the raison d'etre for the theoretical effort undertaken here. 

While current experiments excel in measurements of eigenspectra of compact resonators , 
the periodic orbit theory computations of such bound system spectra are rendered difficult 
by the presence of non-hyperbolic phase space regions. As our primary goal is to derive and 
test rules for replacing wave mechanics by the short-wavelength ray-dynamic trajectories, we 
concentrate here instead on the problem of scattering off cylindrical cavities, for which the 
classical dynamics is fully under control. In the case of one cavity the exact scattering spectrum 
is known [7]. For the multiple-cavities case we generalize the quantum-mechanical (QM) S- 
matrix formalism for A^-disk scattering [8-11], and compute the exact resonances and the 
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Wigner time delays from the full elastodynamic wave-mechanical scattering matrix. We then 
compare the exact results with the corresponding quantities calculated in the short- wavelength 
approximation (SWA), and discover that the QM intuition fails us: the Rayleigh surface waves 
(which have no analog in the QM scattering problem) dominate the low-frequency spectrum 
because of their weights and number, such that already the two-disk elastodynamic scattering 
problem displays chaotic features in this regime in contrast to its QM counter part. 

Elastodynamics. - Consider an infinite slab of an isotropic and homogeneous elastic ma- 
terial {e.g., polyethylene or isotropic quartz) with parallel top and bottom plane boundaries, 
and an in-phase stimulus such that the system behaves quasi-two-dimensionally along the 
slab, with no excitation of or coupling to waves propagating perpendicular to the slab. The 
propagating waves are either the pressure or the shear solutions of the Navier-Cauchy equa- 
tion [12] fi'V^u + (A + /i)V(V • m) + puj^u = 0, where u is a vectorial displacement field, 
A and /i are the Lame constants, p the mass density and lu the frequency. The experiments 
dictate free boundary conditions, with vanishing traction t{u) — where 



t{u) = 



A V 



h. (1) 



Here n is a unit vector normal to the boundary, IL the unit matrix, and T indicates a transpo- 
sition. Elastodynamic waves are vectorial, with the pressure (longitudinal - L) wave propa- 
gating through the bulk with velocity cl = \/ + 2/i)/p, and the shear (transverse - T) wave 
with velocity ct — \/JjJp. Furthermore, the Rayleigh surface waves propagate along plane 
boundaries unattenuated, with the velocity cr determined [12] by the condition 

= (277^ - 1)' - 477^ iv^ - l){v^ - (ct/cl)') with tj = ct/cr. (2) 

When either a shear or pressure plane wave hits a boundary, mode conversion can take 
place [13], with waves of different types emitted at different angles. We now drill one, two, or 
more cylindrical cavities of radius a perpendicularly through the slab. 

One-cavity scattering, the exact spectrum. - Scattering off a single cylindrical cavity is 
separable in angular momentum. The one-cavity S-matrix of elastodynamics, a [2 x 2] matrix 
in {L, T} components, is determined by the free boundary condition for the traction matrices 

SilL = 'Snun' [tiV] • [t(-)] , m = 0, ±1, ±2, • • • . (3) 

For comparison, the corresponding one-disk scattering matrix [11,14] of quantum mechanics 
(QM), for a scalar field with the Dirichlet boundary condition, is expressed in terms of Hankel 
functions, 5'^^/ = — (5„im'H^,^(fca)/Hm ^fca), as function of the wave number k. The traction 
matrices [15] of angular momentum m and superscript Z G {-I-,—, J} result when outgoing 
(-|-), incoming (— ) or regular (J) pressure and shear displacements u in terms of Hankel or 
Bessel functions, Z„i e {H^^^Hm Jm}, arc inserted into the traction 



,2 



(Z)] _ ^ 

(4) 

Index TT e {1,2} labels the 2-d spherical components {r, 0} of the displacement u at the 
cavity while a € {1,2} labels its polarization {L,T} where k^ = ujjca are the corresponding 
wave numbers. The scattering resonances for the one-cavity elastodynamic medium with 
traction-free boundary conditions are determined [7] from Det[t^^] = 0, see eq. 
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Multi-cavity scattering, exact spectrum. - We construct the S-matrix for the iV-cavity 
system foUowing ref. [11]. The determinant of the A^-scatterer S-matrix factorizes, as in QM, 
into the product of the one-cavity determinants and the muhi-scattering contributions [11,15]: 



detSM^MM(^;q;indetSa)^-M. (5) 
^ ' Det[MH] 11 ^ ^ ^ ' 



Here the matrix M = II + A is the inverse of the muhi-scattering matrix. The transfer matrix 
A evolves the displacement of internal angular momentum V = 0, ±1, ±2, • • • and spherical 
component tt' e {1,2} at cavity j' (of radius aj>) to the (/ = 0, ±1,±2, • • • and tt e {1,2}) 
displacement at cavity j (of radius aj) where = 1, • • • , N are the cavity-labels: 



2 2 

a," \ ^ \ ^ 



a," 



= 1 cr' = l 



The matrix [Y^j^'^^^ ]„„i= 5(^(^/H|^;, (fco--Rjj' ) ^xp [iZaj'^ — iZ'(ajj' — tt)] translates the {L, T} 
modes evaluated relative to the origin of cavity j' to the corresponding modes at cavity j. 
Rjji and ajj' are the relative center-to-center distances and angles, respectively [11,15]. The 
multi-scattering resonances of the iV-cavity problem are given by the zeros of Det M(ci;), see 
eq. ^ . For the two-cavity system the two- fold reflection symmetry implies that the determi- 
nant factors into four irreducible representations, with the transfer matrix for each irreducible 
subspace defined on the fundamental domain, a quarter of the full elastodynamic slab [2]. 

Here we present typical numerical results for the fully symmetric Ai subspace of the system 
of two cylindrical cavities. In all presented calculations we take values of the Lame constants 
corresponding to polyethylene [7], and set cl = 1950 m/s, ct = 540 m/s, cr « 513 m/s, and 
take cavities of radius a = 1 cm, center-to-center separation R — 6 cm. The lowest few hundred 
exact Al resonances (hopefully, all A\ resonances in the window [0 < RefcLi < 45,-0.55 < 
ImfcLfl < 0]) determined by the zeros of Det(Il -I- A)|^^ are shown in fig.^ 

One-cavity scattering, ray-dynamic interpretation. - In the one-cavity case a sophis- 
ticated theory of ray dynamics already exists. The main tool is the Sommerfeld- Watson 
transformation which in the QM case replaces a slowly converging partial wave sum by a fast 
converging sum over complex creeping trajectories, first derived by Franz [14]. The pressure 
and shear Franz resonances also exist in elastodynamics [16], but here the one-cavity spectrum 
is dominated by the weakly damped Rayleigh resonances [7], with no QM counterpart. 

In the spirit of Keller's geometrical theory of diffraction [16] we use the one-cavity S-matrix 
to assign a ray-dynamic weight to a segment of the boundary traversed by a Rayleigh 
wave. A circular Rayleigh segment of arc length Ai/ia (with the pertinent angle) has 
the complex- valued weight exp[iA0i^pt(a;)], such that the effective arc length of a Rayleigh 
segment is t^4> v^^iiS) j k'-^ k, t^^aT\, where r\ = ct/cr, see eq. (jSJ. Here we take A(/) times the 
exact wave-mechanical zyR(a;), determined as complex angular momentum by Det[t^^|^^] — 0, 
as the input for the complex-valued action of the Rayleigh trajectory segment. 

Multi-cavity scattering, ray-dynamic interpretation. - In QM the Gutzwiller-Voros Zeta 
function is the semi-classical approximation to Det M , where the connection follows via the 
semi-classical reduction of the traces Tr(A") appearing in the cumulant expansion Det M = 
1+Tr A-[Tr A^ - (Tr A)^] /2-\-- ■ ■ . As shown in refs. [10,11], in the SWA (applicable here from 
/clo ^ 2 onwards) the traces TrA" reduce to the set of ray-dynamic periodic orbits of topologi- 
cal length n, whereas the cumulants become "curvatures" (periodic orbits shadowed by pseudo 



4 



EUROPHYSICS LETTERS 








5 10 15 20 25 30 35 40 45 

Re ki^a 



Fig. 1 - The lowest few hundred exact Ai resonances of the two-cavity elastodynamic scatterer, in the 
complex pressure wave number ki, = lj/cl plane. The labels refer to the various resonance families 
interpreted (together with the lines) in the Discussion section, where the italics specify the orders in 
the cumulant expansion at which the structures are first observed. The second label indicates whether 
the bands are leading (closest to the real axis), additional, subleading, rising, regular or dropping. 



orbits [2]). The shortest geometrical periodic orbits, bouncing between the two cavities, are the 
border orbits of topological length one in the fundamental domain. Their weights, derived from 
the SWA to TtAa,, are of form [11,15] U = -exp[ika{R-2a)]/[y^\K^\{~AoY-^{l- A^'^)] 
with Ao = (i? — a + \J B? — 2Ra)/a, a G {1, 2} = {L, T}. These unstable geometrical orbits 
including their repetitions are summed up in the usual way in terms of the Gutzwiller-Voros 
spectral determinant Det Mgoom. [2]. 

Moreover, Keller's theory yields orbits with unstable geometrical legs and weakly damped 
Rayleigh surface wave arcs circling the cavities. At the topological length one, two (unstable) 
complex Ray leigh- type periodic orbits, an "oval" and "figure eight" contribute [10]. In addition 
to the repeated primary orbits and the "oval-eight" combination, two new (unstable) Rayleigh- 
type orbits contribute to Tr A^, see fig. Eland so on for longer Rayleigh- type orbits, with their 
total number (including repeats) growing exponentially with their topological length. By 
analogy to creeping orbits [17, 18], if gi, i = 1, • • • , n (with qn+i = Qi) are the points along 
a cycle where a Rayleigh arc segment connects to geometric trajectories (which may include 
reflections with/without conversions between pressure and shear rays), the Rayleigh surface 




Fig. 2 - Two new periodic orbits of Rayleigh type of topological length 2 (and "effective" arc radius 
rja > a) in the fundamental domain (and the full space) of the two-cavity system. In the former the 
Rayleigh orbits are further classifiable by the number of contacts with the two symmetry axes [17,18]. 
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Fig. 3 - Wigner time delay for the two-cavity system, the exact Ai resuh versus its ray-dynamic 
expansion to period two, as function of the pressure wave number (times the cavity radius a). 

wave contribution to the spectral determinant is of form 



Here the summation goes over all prime periodic orbits p and their repetition number r. 
[G(gi+i, (7i)]CTi+iCTi (with a^, Oi+i G {L, T}) is the Van-Vleck propagator (including instabilities, 
reflections and mode conversions [15,19]) if qi and Qi+i are connected by a pure geometric 
trajectory, and |[DR(ijj)]cr;_|_icri exp(i (w)) if qi and g^+i are the endpoints of a Rayleigh 
arc segment, with Ui, Ci+i the polarizations of the two attached geometrical legs. The [2 x 2] 
diffraction matrix Dft(a;) is the elastodynamics analogue of the square of the corresponding 
QM diffraction constant [11,17,18]. It is determined by Keller and Karal [16] for the half-plane 
case, and in [15] for the circular cavity. [DR(a;)]CTi_|_iCTi is proportional to exp[—A(a;)CTi+i — A(aj)crJ 
with A((jj)(j = i^Rarccosh[j/R/(aA;o-)] — ('^^o-)^i it is weakly attenuated for the pure 

shear case Ui+i = = T at low frequencies, but strongly attenuated otherwise. 

The "semiclassical" form of the spectral determinant is given by the formal product 
Det Mjgeom. X Det MjRayi., evaluated in the cycle expansion where the classifying topological 
length is equal to the number of geometrical straight legs of the orbits and pseudo orbits [2] . 

The Wigner time delays and cumulant traces are particularly well suited to detailed com- 
parisons of the exact results with the SWA [11]. In fig.|31we plot the exact Wigner "time" delay 
Tci = -^^VcX of the cluster phase shift ??ci(w) — —i\ ln[Det M(w*)'''/Det M(a;)] as function of 
the pressure wave number A^l = ^^/ci^, and compare it to the cycle expansion based on the ge- 
ometrical and complex periodic orbits of topological lengths one and two. The periodic orbit 
approximation is in good agreement with the exact result which can be truncated at second 
cumulant order for the presented A:l values. Finally, in fig. 01 we compare the periods of the 
ray-dynamic periodic orbits with the Fourier peaks of the exact wave-mechanical data by plot- 
ting the moduli of the Fourier transforms of TrA^i^ and Tr for the region 10 < ki^a < 45, 
and the corresponding Fourier transforms of sums over periodic orbits of topological length 
one and two, respectively. The Fourier peaks indeed correspond to the periodic-orbit periods 

= 'iL/cL + ^^t/ct + va^^'in/cT, with li^, li^ and TyaAt/)^^, the (effective) 
geometrical lengths of the pressure, shear and Rayleigh segments, respectively. 
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Fig. 4 - Comparison of the Fourier transforms of the exact traces Tr Aai (left panel) and Tr A^^ 
(right panel) of the two-cavity system with the short-wavelength approximation (SWA) based on 
periodic orbits of topological length one and periodic orbits of topological length two, respectively, in 
the range 10 < khCt < 45. All peaks are identifiable as the periods of the shortest period orbits. 



The SWA captures nearly all qualitative features of the exact calculation. However, the 
1 /ki^a corrections are still not negligible, as can be seen from the small structure close to 0.1 ms 
in fig. 2] (right). This corresponds to a combined LT geometric orbit that vanishes in the SWA, 
as there is no coupling of pressure to shear rays at perpendicular impact. The complex orbits 
visible in fig. 01 have only shear segments coupled to Rayleigh segments, whereas the coupling 
of pressure segments to Rayleigh segments is severely suppressed (since cl ^ cr). 

Discussion. - These observations combined with the cumulant expansion (in detail be- 
yond the scope of this letter) explain the qualitative features of the low-frequency resonance 
spectrum shown in fig. ^ The dotted and dashed lines represent the imaginary parts Im kj^ 
= -iln(Ao)/(i?- 2a) « -0.29/a and Imfcx = -f ln(Ao)cT/cL(i? - 2a) ^ -0.24/a predicted 
by the isolated pressure tj^ and shear orbit, respectively. As the transverse oscillation of 
the tT orbit strongly mixes with the Rayleigh waves of nearly the same wave velocity, the 
irregular band with the leading and subleading resonances arises from the interference of the 
Rayleigh orbits with the ^t orbit. Due to this mixing, there are no resonances on the <t line. 

The Rayleigh orbits include the factor (1 — exp(27rit'R))~^ which arises from the geometrical 
sum of the additional Rayleigh waves around the half-cavity in the fundamental domain. These 
1-cavity structures, which lead to the leading resonances in the 1-cavity case, correspond 
here to poles in detM(fc) which amplify the contributions of the Rayleigh orbits already 
at topological length one, but which also complicate the searches for nearby genuine multi- 
scattering resonances, the zeros of det A/(fc). This is especially the case for the family i-add. 

Below Re/cLO « 32 SWA cycles of at most topological length two are needed to get a 
qualitative fit of the chaotic band labelled as ii-lead, whereas already order one is sufficient for 
the resonances i-lead, before they merge with the rising band ii-ris at about RcfcLO w 4, and 
for the resonances i-add. Below RefcLfl ~ 2, the SWA expansion breaks down and a uniform 
approximation based on a multipole expansion as in Ref. [20] should take over. The dropping 
band, iii-drop, has a QM analog which is generated by Franz' creeping waves. The regular 
family of resonances (see iii-reg and iv-reg) corresponds to the orbit that dominates the QM 
two-disk case [17], but is subdominant here, and that couples neither to the complex Rayleigh 
orbits nor to the shear orbit. When the rising band (iv-ris) crosses the regular band (iii-reg), 
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the cumulant orders three and four are interchanged. It continues as band iii-ris and merges 
with the subleading band ii-sub at about RefcLfl — 35, generating the combined band iii-sub. 

For the fcL window shown in fig-H the periodic orbit sum can be truncated at length four, 
since the spectrum generated from the cumulant expansion to this order does not differ from 
the complete one to the resolution of this figure. This is confirmed by numerical diagonaliza- 
tions of the M x M matrices M^^ (/cl) in the {L, T} x angular momentum space where M has 
to satisfy the bound M > e(cL/cT)|fcL|a {e.g., M > 440 for Re fcL = 45/a) [8,11]. 

Summary and outlook. - We have derived the scattering determinant and calculated, in 
the low-frequency regime, the exact scattering resonances and Wigner time delays for a quasi- 
two-dimensional isotropic and homogeneous elastodynamic slab with two cylindrical cavities. 
Already the physics of this simplest possible multi-scattering system in elastodynamics with 
free boundary conditions is totally different from the one of quantum billiards at low frequen- 
cies: none of the measurable medium excitation {e.g., phase shifts, leading resonances) can be 
understood without the Rayleigh waves which do not have quantum mechanical counterparts. 
For circular boundaries they are barely damped at all [7] and the diffraction constants, which 
link them at the different cavities to the shear waves, are still only weakly attenuated at low 
frequencies. These features are generic for smooth finite-size concave cavities. The pressure 
waves which are the analog of the scalar quantum mechanical ones play only a secondary role. 

A symbolic dynamics in the total space needs to account for cycles patched together 
from 4 kinds of segments: pressure, shear (their sum is the topological length), anti- and 
clockwise Rayleigh, implying an exponentially growing number of interfering periodic orbits. 
Whether the topological increase or the attenuation of the Rayleigh orbits eventually wins 
is still open. Surface orbits of Rayleigh type - with no counterpart in QM - are expected 
in general non-convex elastic resonators. Generalizations to anisotropic media (the highest 
Q-value experiments are performed on single crystals of quartz), and applications of the above 
ray-dynamics techniques to resonator geometries used in experiments remain open problems. 
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